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We present real-time detection measurements of electron tunneling in a graphene quantum dot. 
By counting single electron charging events on the dot, the tunneling process in a graphene con- 
striction and the role of localized states are studied in detail. In the regime of low charge detector 
bias we see only a single time-dependent process in the tunneling rate which can be modeled using 
a Fermi-broadened energy distribution of the carriers in the lead. We find a non-monotonic gate 
dependence of the tunneling coupling attributed to the formation of localized states in the constric- 
tion. Increasing the detector bias above Vb = 2 mV results in an increase of the dot-lead transition 
rate related to back-action of the charge detector current on the dot. 

PACS numbers: 72.80.Vp, 73.63.Kv, 73.50.Td, 73.23.Hk 



I. INTRODUCTION 

The high sensitivity of a quantum point contact or 
a single electron transistor to its electrostatic environ- 
ment is widely used as a powerful tool to detect the 
electron occupation in semiconductor-based quantum dot 
structures.^ ^ Time resolved charge detection,^ offers 
the possibility to measure extremely small currents by 
counting single electron transitions^ and enables e.g. 
the extraction of detailed electron tunneling statistics 
and probing of electron-electron correlations.^ Time- 
resolution further allows for single shot read out of spin- 
qubits after spin to charge conversion with potential ap- 
plications in future quantum information processors.^^ 

With the rise of two-dimensional graphene^^ a fas- 
cinating new mesoscopic material for transport exper- 
iments has become available with the promise of long 
spin coherence times. First nanostructures have been 
made by etc hing graphene into narrow constrictions 
(nanoribbons).I^^I^In these devices transport around the 
charge neutrality point is suppressed. Short constric- 
tions have been successfully used as barriers for graphene 
quantum dots.l^^K^ Several experiments have been per- 
formed on graphene quantum dots including the obser- 
vation of excited states in single^^ and double quan- 
tum dots^^ and the investigation of orbital^^ and spin 
state4^ around the electron- hole crossover. An addi- 
tional nanoribbon placed in close proximity to the dot 
can be used to detect the number of electrons on the dot^^ 
(see also Ref. [27 1) . In contrast to charge detection with 
a quantum point contact (QPC) tuned below the lowest 
conductance plateau, highly charge-sensitive resonances 
in the detector nanoribbon are used as a sensor, similar 
to charge detection with a single electron transistor. [2^1^ 
These resonances arise from localization of charge car- 
riers due to strong potential fluctuations in the disor- 
dered nanoribbon. 29 33 Yiere we investigate tunneling 
through a graphene quantum dot lead in a time-resolved 
way by counting individual charging events with such an 
integrated graphene charge detector. The time resolution 
allows for a deepened analysis of the tunneling properties 
of a graphene constriction compared to the time-averaged 



casfPia. 

The paper is organized as follows: we first briefly char- 
acterize the charge detector and the quantum dot in 
Sec. II. In Sec. Ill time-resolved charge detection is pre- 
sented. It is shown that despite inherent resonances and 
transmission modulation in graphene barriers, transport 
can be well understood within a conventional model. For 
a particular gate-voltage regime an asymmetric double 
barrier model is successfully used to describe the trans- 
mission through a graphene constriction. While increas- 
ing the bias, back-action of the detector on the dot is 
observed. In Sec. IV the performance of the detector is 
discussed. A short summary of the results is given in 
Sec. V. 



II. DEVICE CHARACTERIZATION 

The structure is carved out of a mechanically exfoliated 
graphene flake by reactive ion etching in an argon-oxygen 
plasma. Details of the fabrication process can be found 
in Ref. 34, The all-graphene sample [shown in Fig. 1(a)] 
consists of a 95x70 nm quantum dot (qd) with three lat- 
eral graphene gates: left gate (Ig), plunger gate (pg) and 
right gate (rg) [see Fig. 1(a)]. A nanoribbon is used as 
a charge detector (cd) with one additional lateral charge 
detector gate (cdg) for tuning the detector to the regime 
of highest sensitivity. The highly p-doped Si substrate, 
isolated by 295 nm Si02, is used as a back gate to tune 
the Fermi energy. The constrictions connecting the dot 
to the leads are only 15 nm and 20 nm wide to achieve 
very low tunneling rates of the order of 1 Hz- 100 Hz. The 
narrower constriction turned out to be completely isolat- 
ing and no carrier tunneling between dot and source was 
observed. All measurements were performed in a ^He 
cryostat at a base temperature of T = 1.7 K. 

By changing the back-gate voltage the charge detector 
can be tuned from hole dominated transport at Vbg < 
— 11.5 V to electron dominated transport at Vbg > 0.5 V 
[see Fig. 1(b)]. In between the conductance is pinched 
off and shows resonances which are the signatures for 
a "transport gap "^^^-^^ where the conductance is gov- 
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FIG. 1: (Color online) (a) SFM micrograph of of the graphene 
quantum dot with source (s), drain (d) and lateral gates (pg, 
Ig, rg) as well as the charge detector (cd) and its gate (cdg) 
(b) Measurement of the conductance through the CD as a 
function of back gate voltage Vbg at a bias of V^^ — 10 mV. 
The cones indicate hole (left) and electron (right) transport, 
whereas the reddish area marks the transport gap. (c) Single 
resonance in the charge detector conductance recorded as a 
function of at V^^f '^^ = -5.66 V with V^"^ = 7mV. The 
effective Vbg contains the contribution from the back gate 
Vbg = —2.083 V and the influence of the charge detector gate 
Vcdg — —11.56 V (as in the other measurements shown in 
the paper while all other gates including V^^^ are at V ). 
The abrupt changes in Gcd arise from charging the dot with 
additional electrons. 



erned by localization and Coulomb blockade induced by 
disorder The center of the gap at Fbg ~ 5 V is offset 
from Vbg = indicating an overall negative doping of the 
charge detector. 

Fig. 1(c) shows a close-up of a resonance recorded by 
changing the left gate voltage V\g while all other gate 
voltages are fixed in the gate regime indicated by the 
dashed line (see caption). Due to the high sensitivity 
of the detector in this gate configuration, this regime is 
also used in the following measurements. The kinks on 
the resonance originate from charging the dot electron by 
electron as V\g is increased. Whenever an additional elec- 
tron is loaded to the dot the electrostatic potential in the 
detector changes abruptly, which results in an abrupt in- 
crease (reduction) of the detector conductance if the po- 
tential energy is above (below) the resonance condition in 
the detector. The gate voltage shift /SV\^ ^ 0.13 V of the 
curve at each step corresponds to the potential change in 
the nanoribbon due to the additional dot electron. 

Note that the spatial locations of the charging events 
are identified by analyzing the influence of the gates 
on the jumps and resonances in Gcd (not shown here, 
see Ref. 26). The influence is characterized by a lever 
arm ag which relates the gate voltage Vg to the in- 
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TABLE I: Lever arms of the different gates on abrupt quan- 
tum dot conductance changes and the resonances in the 
charge detector. The lever arms are extracted by compar- 
ing the influence of each gate relative to the back gate. The 
absolute plunger gate lever arm on the (dot-) charging events 
is extracted from dot-bias dependent measurements by as- 
suming symmetric coupling to source and drain. The other 
absolute dot lever- arms are calculated based on this estimate. 

duced potential change in device x. The result is tab- 
ulated in Tab. 1 and shows a similar influence of the 
right and left gate on the abrupt changes of Gcd ^ind 
superior coupling to the plunger gate. Such a behav- 
ior is expected if the abrupt changes result from charg- 
ing the dot. Moreover, the absolute values of the lever 
arms can be estimated using dot bias spectroscopy by 
assuming symmetric electrostatic coupling to source and 
drain, revealing an addition energy of £^c = 19 meV (not 
shown) . This is in good agreement with the charging en- 
ergy £;^isk _ ^ ji^/^^^^^ ~ 20 meV for an isolated disk of 
diameter d = 90 nm, expected (as an upper limit) from 
the dimensions of the island. The spacing of the jumps in 
gate voltage are then related to an absolute plunger gate 
lever arm ol^ = 0.10, based on which the other absolute 
lever arms in Tab. 1 are estimated. 

The height of the step A/gtep and hence the magni- 
tude of the detector signal depends (i) on the induced 
change of the detector potential {ea^^AV\g ~ 1 meV) 
given by the electrostatic coupling between the dot and 
the detector and (ii) on the sensitivity of the detector 
conductance to potential changes (steepness of the de- 
tector current measured as a function of gate voltage). 
Here, afg is estimated from geometrical considerations 

to be a^^^ - a^^g ^ 0.02. 



III. TIME RESOLVED CHARGE DETECTION 

A. Time dependent detector current around a step 

In the following we analyze the detector conductance 
around such a step in real time. Fig. 2(a) shows a zoom 
of the time averaged current across a step measured with 
a typical integration time of 0.2 s. Here a lower bias of 
V^^ = 0.5 mV is applied. We don't observe a single step 
(see gray line) but the signal is noisy and the transition 
is smeared out. These are indications that the timescale 
of the transition is comparable to the measurement in- 
tegration time. Indeed, by measuring a time trace at 
position (b) [see Fig. 2(b)] we observe a two-level ran- 
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FIG. 2: (Color online), (a) Time averaged current through 
the charge detector as a function of Mg while scanning over 
a similar resonance in the QD as shown in Fig. 1(c). The 
step indicates a single change from N to N+1 electrons on 
the dot at V^^ = 0.5 mV. (b) Time resolved current through 
the charge detector taken at the point labeled with (b) in (a). 
The time an electron needs to tunnel into or out of the dot is 
marked with Tin and Tout respectively, (c) Histogram of the 
led values for the whole time trace of 60 s. (d-g) Time traces 
taken at different V^g marked in (a), show a gradual change of 
the dot electron number from N (lower level) to N+1 electrons 
(upper level). 



dom signal switching in intervals of around 0.5 s. The 
two-level signal shows a large signal to noise ratio (SNR) 
of A/step/ (^noise) ~ 30 at a measurement system band- 
width of ^ 400 Hz. The noise (/noise) is defined as the 
variance Y^Var(/) of the detector current on each of the 
two current levels. The two levels indicate whether there 
are N electrons (lower) or + 1 electrons (higher level) 
on the dot and hence allow real time detection of sin- 
gle charge carriers tunneling on and off the dot. The 
corresponding dwell times Tin and Tout are indicated in 
the figure. The histogram in Fig. 2(c) shows the distri- 
bution of the current for a 60 s time trace with a clear 
separation of the two states. The slight asymmetry in 
the distribution of the upper level is attributed to an ad- 
ditional weakly coupled charge fluctuator present in the 
+ 1 state. 

By tuning the dot potential away from the resonance 
condition the occupation probability of the lower N [up- 
per N level is reduced as seen in Fig. 2(e,d) [f,g]. In 
Fig. 2(g) for example the chemical potential for the N-\-l 
transition lies below the chemical potential in the lead. 
Therefore the dwell time for the empty state Tin is much 
smaller than Tout- 
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FIG. 3: (Color online), (a) Distribution of the dwell times 
recorded at a position where the dot chemical potential is 
slightly above the drain chemical potential [see schematic and 
indicated by the arrow in (c)]. The distribution shows an ex- 
ponential behavior described in Eq. ^ with Fin = 47 Hz and 
Tout = 102 Hz. (b) Occupation probability for the N+1 elec- 
trons dot state P(N +1) as a function of Vig. The solid 
line is a fit to 1 — f{E) with f the Fermi function. Assum- 
ing T = 1.7 K a Vig lever arm of aig = 0.077e is obtained. 
For comparison the dashed line shows the result obtained 
in (c). Here the events per second are plotted for varying 
V\g. Comparing these event rates with Eq. Q the tunnel- 
ing coupling and the lever arm are given as F = 143 Hz 
and aig = 0.067e (solid line). The difference of the lever 
arms extracted from (b) and (c) (dashed line) is attributed 
to the limited statistics, (d) The data from (c) is fitted in 
addition with a tunnel-coupling broadened Lorentzian line 
shape oc F/ (^aig(^Vig)^ + F^) (dashed line) and plotted in 
logarithmic scale, confirming the thermal broadening of the 
resonance. This measurement has been conducted at a dot 
resonance close to the one analyzed in Fig. 2 with faster tun- 
neling rates [see Fig. 4(a)]. 



B. Quantitative analysis of time traces 

The following analysis follows closely recent work on 
similar systems realized in semiconductors.^* ^^ * ^^ * A quan- 
titative analysis of time traces reveals information about 
the number of dot-levels participating in transport, the 
tunneling coupling, the carrier temperature and distri- 
bution in the leads and the individual tunneling rates. 
For this analysis the barrier is slightly opened to obtain 
larger count rates and therefore improved statistics. In 
Fig. 3(a) the distribution of the dwell times is plotted 
for a situation where the chemical potential of the dot 
state /iAT+i is slightly above the chemical potential of the 
drain /id (see schematic). The dwell times are exponen- 
tially distributed indicating that only a single dot level 
is involved in the transport process. Hence the probabil- 
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ity density Pin/out (^) which is the number of counts with 
dweh time Tin/out = t normahzed by the total number of 
counts is given by 



Pin/ out {t)dt = 



1 



('^in/out) 



exp 



t 



('^in/out) 



\dt. (1) 



From a fit of this equation to the data in Fig. 3(a) tunnel- 
ing rates Fin = = 47 Hz and Fout = — 
are extracted. 

As shown in Fig. 2(d-g) these rates change by tuning 
the potential of the dot. Fig. 3(b) shows the occupa- 
tion probability of the N+1 state changing from to 1 
while lowering the chemical potential of the dot for in- 
creasing V\g (see schematics). If we assume gate voltage- 
independent tunneling coupling to the lead the occupa- 
tion probability is determined by the distribution of the 
charge carriers in the lead. 

We assume a Fermi distribution in the lead / = 
[1 + exp (A/i//cTe)]~^ with Aja the difference between 
the chemical potential of the dot and the lead. By tun- 
ing the left gate voltage A/i is changed according to 
A/i = eaig{Vig — Vres) where aig is the lever arm of 
the left gate on the dot (the small lever arm on the 
lead is neglected). The occupation probability is then 
given as P/v+i(A/i) = 1 — /(A/i) and we can extract 
a\g ~ 0.077 ± 0.01 under the assumption that the elec- 
tron temperature Tq is equal to the bath temperature of 
1.7 K. The uncertainty arises from the limited statistics 
around the crossover point and is reduced by analyzing 
the number of tunneling events (see below). The lever 
arm is comparable to a\g ~ 0.06 estimated from the dot 
bias dependence in Sec. IL 

Combining single-level transport and constant tunnel 
coupling it is possible to extract the tunnel coupling F by 
counting the number of tunneling events [see Fig. 3(c)]. 
The event rate Tq for tunneling-in is given by 



1 



(Tin) + (Tout) 



r-/(i-/). 



(2) 



The best fit to the data yields F = 143 Hz and aig ^ 
0.067 ±0.005. The fit obtained in (b) is plotted for com- 
parison as a dashed line in (c) and vice versa in (b). 
The discrepancy between the two fits is explained by the 
larger influence of the data points around the resonance 
in (b) compared to (c). 

In Fig. 3(d), a Lorentzian broadened line shape 
(dashed) and the thermally broadened fit (solid) of the 
number of events are plotted on a logarithmic scale for 
comparison. As expected for low tunneling coupling the 
thermally broadened line shape describes the data much 
more accurately. 



C. Changing the tunneling barrier 

By changing the corresponding s ide ga te we can tune 
the barrier and the tunnel coupling .^^^ Fig. 4(a) shows 
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FIG. 4: (Color online), (a) Number of events per second 
measured as a function of Vig and Kg while crossing a dot 
resonance. A single trace taken at Kg = —0.37 V has been 
analyzed in Fig. 3. (b) Tunneling coupling as a function of 
Kg (and ^ —Kg) plotted in logarithmic scale. The red circles 
are obtained from the data shown in (a) (regime A), whereas 
the blue squares are deduced from a similar measurement in 
a slightly shifted gate regime B around Kg = V and Kg = 
1.75 V after a small charge rearrangement (lOx magnified for 
clarity). The non-monotonic behavior can be explained by 
modeling the constriction as asymmetric double barrier with 
lorentzian shaped resonances from the weak barrier (/iFa = 
6.2 meV and HTb =6.4 meV) and an exponential suppression 
due to the strong barrier (see text for details), (c) Peak width 
plotted as inverse lever arm 1 / arg and FWHM (right scale) as 
a function of Kg in both regimes with an average arg = 0.062. 
Here V^^ = 0.5 mV and the counting time is 60 s per point. 



the number of events per second as a function of Kg and 
Kg- The measurement analyzed in Fig. 3 is a cut at 
Kg = —0.37 V. The diagonal line corresponds to the res- 
onance condition where an additional carrier is loaded 
onto the dot (N N+1). The potential of the dot is 
affected equally by both gates as expected from the ge- 
ometry [Fig. 1(a)]. In addition a change in the number 
of counts is observed especially at lower Kg- This change 
in the barrier transmission F can be seen more easily in 
Fig. 4(b) where the corresponding tunneling coupling F 
[obtained by fitting to Eq. ([2|] is plotted as a function of 
Kg (A, red circles). The same behaviour is also observed 
in a second measurement (B, blue squares) taken in the 
same gate regime after a small charge rearrangment and 
offset with a factor of 10 for clarity. While the tunneling 
coupling varies strongly by changing Kg the full-width- 
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half- maximum (FWHM), which is inversely proportional 
to the lever arm is approximately constant with an aver- 
age aig = 0.062 at T = 1.7 K [see Fig. 4(c)]. This aver- 
aged lever arm is in good agreement with the aig ^0.06 
obtained by varying the dot bias. 

Unlike GaAs based quantum dots where tunneling 
rates tend to increase monotonically with gating due to 
depletion of the electron gas^^, the non-monotonic ex- 
ponential changes of the tunneling rate in our graphene 
nanostructure is an indication for the presence of reso- 
nances in the constrictions.^^ It is important to note 
that we see no charging of constriction resonances in this 
measurement. However, it is possible to tune the device 
into a regime where the typical hexagon pattern of a dou- 
ble dot is measured while changing versus V\^^ and 
signatures of a second localized state in the barrier can 
be observed. If charging of well localized parasitic reso- 
nances is slow enough even additional small steps in the 
counting signal are observable in those regimes. 

In the gate regime investigated in this paper, the influ- 
ence of localized states on the dot energy is negligible but 
still the barrier transmission is modulated. Such a behav- 
ior might occur if the additional localized state is strongly 
coupled to the lead but only weakly to the dot. This sit- 
uation can be modeled with a one-dimensional asymmet- 
ric double barrier with tunneling coupling F^^ F/,. For 
noninteracting electrons in the case of kT <C /iF <C A 
(with A the level spacing) the total transmission is given 
aPMH 



'(mV) 
4 



tot 



(rR/2)2 + (Ef - EpY 



(3) 



if we assume Fl,r to be independent of the resonance 
p. In this limit the Lorentzian shape is caused by the 
weak barrier with strong coupling while the overall am- 
plitude is determined by the weak tunneling coupling of 
the strong barrier. For the strong barrier we assume an 
exponential dependence of Fl on the gate voltage while 
the gate dependence of the weak barrier is neglected. The 
dashed lines in Fig. 4(b) are the corresponding fits with 
two resonances {p = 1,2). The extracted tunneling cou- 
pling of the weak barrier is similarly strong in both mea- 
surements with Fr A = 6.2 meV and Fr b = 6.4 meV 
assuming a typical lever arm of arg,ioc = 0-1 on the 
localized state. For the left barrier the WKB result 
with a linearized exponential Fl = Fq exp [—hz{dU — SE)] 
is used.^ The details of the barrier are described by 
Fq and while SE {SU) describes small perturbations 
of the dot energy (barrier potential). By keeping the 
dot energy constant and assuming a linear gate depen- 
dence of the barrier potential we get Fl = Fq exp [— /3Kg] 
where /3 = K.a^^^^^^^. For the two measurements we ob- 
tain Fo,A = 0.1 Hz, Pa = Q and Fo,b = 0.2 Hz, 
= 11 V~^. The different parameters in the two 
regimes are attributed to a change of the left barrier po- 
tential between the two measurements. 
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FIG. 5: (Color online), (a) Signal to noise ratio for increasing 
charge detector bias V^^. Every data point is obtained from 
the average step height (signal) and the average noise on the 
two levels from 60 time traces (Vig = 1.735 — 1.79 V) each 
25 seconds long at Vrg = 0. (b) Signal and noise versus 
plotted independently. The bold arrow marks the onset of 
the stronger increase in noise giving rise for the saturation of 
the signal to noise ratio in (a). (c,d) Dependence of the dot 
events on the detector bias. Here the measurement B shown in 
Fig. 4(b,c) for ^ = 0.5 mV is repeated for different values of 
V^"^ with -150 mV< Vrg < 25 mV. In (c) the relative change 
of the tunneling rate F is plotted. In order to compensate for 
variations in the tunneling rate by changing Kg , Fq is defined 
as the average tunneling rate for V^^ < 2 mV for each Vrg. In 
(d) the FWHM of the peaks [as shown in Fig. 3(c)] averaged 
over the different Vrg is plotted for increasing V^^. 



D. Detector bias dependence and back-action 

In the measurements shown so far the bias in the 
charge detector was kept low {V^^ = 500 /iV) in order 
to prevent back-action of the detector on the dot.^^^^' 
On the other hand a higher charge detector bias leads to 
an increase of the signal (-step). In order to maximize 
the performance, the SNR is investigated as a function 
of charge-detector bias in Fig. 5(a,b). The SNR is maxi- 
mized for 1^^^ = 2 mV and gets smaller for V-^^ > 2 mV 
due to a stronger increase of the noise [bold arrow in 
Fig. 5(b)], as can be seen in Fig. 5(b), where the signal 
and the noise are separately plotted. This higher noise 
is correlated with an increase in tunneling events and a 
broadening of the dot event peak [see Fig. 5(c,d)]. The 
data is obtained by recording a left-gate right-gate map 
of the charging events [such as Fig. 4(a)] for different V^^. 
The extracted FWHM of the peak and tunneling rates are 
averaged over 15 right gate values (Kg = [—150 — 25 mV]) 
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in regime B (see squares in Fig. 4(b,c) where the tunnel- 
ing couphng and the FWHM are shown as a function of 
Vrg at V^^ = 0.5 mV). In order to account for variations 
of the tunnehng rate with Vi-g, an average tunnehng rate 
To (Kg) for V^^ < 2 mV is defined for each Vrg value. 
In Fig. 5(c) the average rate F/Fq increases up to 40% 
from V^^ < 2 mV to V^^ = 7 mV. Note also the increase 
of the standard deviation of the average value, reflected 
in the errorbar. In Fig. 5(d) the FWHM is calculated 
from the peak width using the leverarm aig = 0.06 ob- 
tained from Fig. 4(c). Similar to the tunneling rate the 
peak width depends approximately linearly on the detec- 
tor bias above V^^ = 2 mV. 

Due to the correlation with the noise in the charge de- 
tector, we attribute the back-action from the detector on 
the dot to arise main ly fro m shot noise generated in the 
detector constrict ion. ^^^ESI Photon emission and absorp- 
tion is rather easy in graphene due to the linear, zero- 
bandga p elec tronic dispersion. Heating due to acoustic 
phonon^^^^^is less plausible because the phonons have to 
couple via the Si02 substrate over a different material^. 
In addition, graphene has a very high thermal conductiv- 
ity 5000 W/m/K)49 compared to SiOs (1.3 W/K/m, 
both measured at T = 300 K) and therefore the gen- 
erated heat is expected to thermalize in the leads of the 
graphene constriction rather than to heat the dot lead via 
the oxide. However, in order to clarify this mechanism 
further experiments including double quantum dots with 
frequency resolved absorption measurements^^ or more 
involved studies of the detector-current dependence at 
constant detector bias are desirable. 



IV. DETECTOR PERFORMANCE 

In the presented device a current step of up to 2 nA 
with a relevant noise spectrum below 2 pA/^/Hz is mea- 
sured at a detector bandwidth of 4 kHz {V^^ = 1.7 mV) 
(see Fig. 6). This corresponds to a charge sensitivity of 
the detector better than 10~^e/V^Hz comparable to what 
has been reported in Ref. 51 for a GaAs QPC detector. 

Although in principle a bandwidth of up to 30 kHz 
can be achieved with a conventional room temperature 
amplifler setup, as has been shown with GaAs QPC 
detectors j ^ l ^ l ^^ l the bandwidth is limited in the presented 
measurements to below 1 kHz to ensure sufficient SNR 
for counting. 

The first limitation is posed by the noise of the sys- 
tem. Usually the noise spectrum of a charge detector in 
such a setup is dominated by intrinsic 1 / f- noise from the 
sample at frequencies / < 1 kHz and by amplifier noise 
at higher frequencies. In our sample the 1/f-noise of the 
device is rough ly thr ee times larger compared with GaAs 
QPC detectors.l^^The negative infiuence on the SNR is 
limited by reducing the bandwidth with a larger feedback 
resistance. Concerning the amplifier noise at higher fre- 
quencies, the comparably high resistance of our detector 
{Red = 500 kQ compared to i^^d^ ~ QPC-based 
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FIG. 6: (Color online), (a) Current- and charge-noise versus 
frequency for different feedback resistors = 10 MQ (blue, 
bandwidth 4 kHz) and 100 MQ (red, bandwidth 400 Hz) mea- 
sured at a bias V^"^ = 1.7 mV with optimal signal-to- noise 
ratio {Red = 570 kQ). The noise is obtained by averaging the 
noise spectra of 10 time-traces, each 0.5 s long (containing 
no dot charging events). The total modeled noise (black) is 
mainly determined by 1/f-noise of the sample (long dashes) 
and the amplifier noise (short dashes). The thermal noise of 
the feedback resistor is lower in magnitude (dash-dotted). On 
the right scale the charge noise is plotted. It is obtained by di- 
viding the noise by the step height of the signal /step = 1.9 nA 
induced by charging the dot with an additional electron. The 
value below 10~^e/ VHz above / = 1 kHz is comparable with 
the resolution shown in Ref. 1511 



detectorJ^^I^ leads to a weaker signal amplification and 
hence the amplifier noise gets proportionally more im- 
portant. Here an increase of the feedback resistance is 
beneficial as well, as the amplification is restored. 

A second limitation for the SNR is related to the gen- 
eral difficulty using SET-based charge detectors to main- 
tain the working point with optimal sensitivity. In prin- 
ciple it is possible to compensate the infiuence of any 
gate on the charge detector with the charge-detector gate. 
However, the change of Vcdg induced additional charge 
fiuctuations in the detector and was therefore kept at a 
constant value during measurements. 

A more general issue for time resolved charge detection 
in graphene is the limited tunability of the dot barriers. 
Although the current through etched graphene quantum 
dots can be tuned over several orders of magnitude the 
tunability of the tunneling barrier is limited by the width 
of the constriction. 

Despite the mentioned issues with the current imple- 
mentation, we could show high sensitivity with signal 
changes of more than 50% at moderate cryogenic temper- 
atures of T = 1.7 K. This is because graphene offers the 
possibility for very strong electrostatic coupling between 
charge detector and quantum dot, because the spacing 
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between them can in principle be made even smaller than 
the 30 nm used in this device. The coupling could be 
even further increased by making use of the monoatomic 
thickness of graphene and vertically stacking dot and 
detector. In addition, the bandwidth can be improved 
using a low te mpe rature amplifie J ^^ * ^^* and / or a radio fre- 
quency setup j^l^^H^ where a bandwidth of 1 MHz has 
been shown with a charge sensitivity of < 2 • 10~^e/ ^/Rz 
in Ref . E 

V. SUMMARY 

We demonstrated time-resolved charge carrier detec- 
tion in a graphene quantum dot with high charge sen- 
sitivity of the detector due to its close proximity to the 
quantum dot. For the measurements recorded at low 
detector bias, the tunneling rate can be modeled con- 
ventionally by a single time-dependent process with a 



temperature broadened energy distribution of carriers in 
the lead. Gating of the tunneling barrier reveals a non- 
monotonic gate dependence of the tunneling coupling by 
counting individual charging events. This behavior is at- 
tributed to resonant tunneling through localized states in 
the barrier strongly coupled to the lead with hV ^ 6 meV. 
For detector bias V^^ > 2 mV we see a symmetric broad- 
ening of the energy distribution of the tunneling events 
accompanied by an increase in the detector noise. This 
back-action is attributed to shot noise in the charge de- 
tector, as graphene offers a high thermal conductivity in 
contrast to Si02 together with the ease of photon emis- 
sion and absorption in an a priori zero-bandgap material. 
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